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ABSTRACT 

We have carried out a pilot study for the SCUBA-2 'All-Sky' Survey, SASSy, a wide and 
shallow mapping project at 850 /im, designed to find rare objects, both Galactic and 
extragalactic. Two distinct sets of exploratory observations were undertaken, and used 
to test the SASSy approach and data reduction pipeline. The first was a 0.5° x 0.5° map 
around the nearby galaxy NGC 2559. The galaxy was easily detected at 156 mJy, but 
no other convincing sources are present in the map. Comparison with other galaxies 
with similar wavelength coverage indicates that NGC 2559 has relatively warm dust. 
The second observations cover 1 deg 2 around the W5-E H II region. As well as diffuse 
structure in the map, a filtering approach was able to extract 27 compact sources with 
signal-to- noise greater than 6. By matching with data at other wavelengths we can see 
that the SCUBA-2 data can be used to discriminate the colder cores. Together these 
observations show that the SASSy project will be able to meet its original goals of 
detecting new bright sources which will be ideal for follow-up observations with other 
facilities. 

Key words: surveys - submillimetre: galaxies - submillimetre: stars 



1 INTRODUCTION 

The millimetre and sub-millimetre parts of the electromag- 
netic spectrum directly probe the cold Universe. The sub- 
mm window specifically allows us to study the youngest 
phases of star formation in our Galaxy, and the dusty, most 
prodigiously star- forming galaxies at high redshift. Despite 
this strong motivation, the sub-mm sky still remains poorly 
surveyed. The SCUBA-2 'All-Sky' Survey, or SASSjQ, is 
a James Clerk Maxwell Telescope (JCMT) Legacy Survey 
(JLS) project designed to redress this balance and exploit 
the rapid mapping capability of SCUBA-2 to ultimately map 
a large portion of the sky visible from the JCMT to an an- 
gular resolution of 14 arcsec at 850 /im. The target point 
source rms level is 30 mJy. 

The benefits of such a wide-field survey are many, rang- 
ing from a complete census of infrared dark clouds (IRDCs) 
to the potential discovery of some of the most luminous high - 
redshift galaxies in the Universe (|Thompson et al l 2007). 
The approved phase of SASSy consists of two distinct parts: 
a strip covering the Galactic Plane which is visible from 
Hawaii; and a 'Pole-to-Pole' strip perpendicular to this and 
designed to pass through the Galactic and Ecliptic North 
Poles. These observations will be carried out in 'Grade 4' 
weather conditions (0.12 < T225 < 0.2, where T225 is the sky 
opacity at 225 GHz, as measured by the CSO radiometer), 
i.e. essentially when the atmosphere is too opaque to enable 
useful observations of fainter objects with SCUBA-2. The 
450 fim data are therefore expected to be of marginal value, 
and the survey is entirely designed to make large maps at 
850 (jm, with a target sensitivity of 30 mJy, achieved using 
a fast scanning speed. 

SASSy will be able to build on the success of IRAS at 
one decade shorter wavelengths. It is also complementary to 
several more recent wide surveys, namely WISE and Akari, 
at near- and mid-IR wavelengths, and surveys with the Her- 
schel and Planck satellites in the submillimetre (hereafter 
submm). Early results from Herschel have already demon- 
strated that bright lensed galaxies can be selected at submm 
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1 Alternatively known as the SCUBA-2 Ambitious Sky Survey. 



wavelengths iNegrello fc et al.l l|2010l ). and also that poten- 
tially protostellar cores can als o be found using Herschel 
IWard- Thompson fc et all |2010l ). The Planck Early Release 
Compact Source Catalogue will contain an all-sky list of 
sources which includes the 850 fim channel. However, the 
resolution of SASSy will be 20 times better. Hence it is clear 
that if SCUBA-2 is able to map rapidly enough while main- 
taining its nominal beamsize and achieving the required sen- 
sitivity, then SASSy will be able to meet its science goals. As 
we show below, these preliminary SCUBA-2 data lend con- 
fidence that these requirements can be met and that SASSy 
will discover many new and interesting sources. 

In the next section we describe the two sets of obser- 
vations which were made as a pilot study for SASSy, one 
Galactic field and one extragalactic field. In Section [3] we 
describe how we reduced the data and extracted sources. 
We then discuss the properties of the sources that we found 
and end with some conclusions about SASSy in general. 



2 OBSERVATIONS 

SCUBA-2 (|Holland et alj|2003l , [200l ) is the successor to the 
Submillimetre Com mon User Bolometer Array (SCUBA, 
iHolland et al.l [l999), which operated successfully on the 
JCMT from 1997-2005. SCUBA-2 has been designed to 
be hundreds of time faster than SCUBA for mapping the 
sky in the same two primary bands, 450 /im and 850 /im. 
The observatory offered a period of 'Shared Risks Observ- 
ing' programme (hereafter S2SRO), in which relatively short 
programmes were carried out with a partially-commissioned 
version of the instrument, having one (of four) sub-array 
available at each of the two bands. Unless stated otherwise, 
all maps and detections in this paper are at 850 /J,m, since 
this is the primary wavelength of interest for SASSy. 

As part of S2SRO there were two separate sets of obser- 
vations carried out for SASSy, which we now describe. The 
SCUBA-2 array of 32 x 40 transition-edge sensitive bolome- 
ters has a footprint of about 3 arcmin square on the sky (it 
will double in linear size when the focal plane is fully popu- 
lated with sub-arrays) on the sky, and takes data at a sam- 
pling rate of approximately 200 Hz while scanning across the 
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Figure 1. Smoothed signal-to-noise map of the 'extragalactic' 
field, being about 0.5° across. NGC 2559 is located at the centre of 
this map, which is dominated by artefacts of roughly the SCUBA- 
2 array size. 
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Figure 2. Matched- filtered signal-to- noise map of the 'extra- 
galactic' field. NGC 2559 is located at the centre of this map, 
and is clearly detected. Other peaks in this map (near the north- 
east and south-east corners) appear to be just noise excursions. 



2.1 Extragalactic observations 

The first set of observations planned was imaging of a region 
measuring 0.5° x 2°, with one end centred on NGC 2559, a 
nearby galaxy, and the other end extending into the Galactic 
Plane cirrus. NGC 2559 is a spiral galaxy at a distance of 
20.8 Mpc (from the systemic velocity and assuming a Hubble 
constant of 75 kms" 1 ), with morphological type SB(s)bc pec 
|de Vaucouleurs et al.|[l99l} ). It was chosen as a target in this 
study since it is the IR-brightest nearby galaxy that lies in 
the SASSy area, and had not been previously observed by 
SCUBA. This section of sky was broken up into four 0.5° x 
0.5° tiles, of which only two were observed during the S2SRO 
period - these were the outer ends of the strip, and hence 
are not contiguous. We refer to this as the 'extragalactic' 
pilot survey, since (although relatively close to the Galactic 
Plane) the field was chosen to contain a bright IRAS galaxy 
which had not previously been observed at submillimetre 
wavelengths. 

The region containing NGC 2559 was observed on 2010 
February 27 and March 14, for a total of 167 minutes of 
observing yielding an average integration time of 27 seconds 
per pixel. The average optical depths at 225 GHz for the 
two nights were T225 = 0.097 and 0.154, the average noise- 
equivalent flux densities were 160 and 260 mjys" 1/2 , and 
the telescope scan rates were 240 and 360 arcsecs -1 , re- 
spectively. Data were reduced using the SMURF map-maker 
(Sub-Millimetre Us er Reduction Facility. TChapin et al.ll20l61 : 
Ijenness et alJl2O10L Chapin et al. in preparation), which we 
describe in more detail in the next section. The resulting 
raw map has a noise of 38mJybeam _1 determined using 
the produced noise map. 

The region about 1° away, containing known Galactic 
cirrus, was observed for a total of 105 minutes, resulting in 
an average integration time of 16 seconds p er pixel. Both 
'pong ' and 'rotating pong' scan strategies (jKacklev et al.l 
2010B were used for the NGC 2559 field and only the rotat- 
ing pong strategy for the cirrus field. Although NGC 2559 
is readily detected in the resulting map of this area, we are 
unable to identify any Galactic cirrus within the other map. 
This is due to a combination of the relatively high noise 
level, and the difficulty in detecting extended diffuse struc- 
ture (largely removed in the reduction process). Therefore 
for the remainder of this paper, of the two extragalactic ob- 
servations, we restrict our attention to the map containing 
the galaxy NGC 2559. 



2.2 Galactic Observations 

The second distinct set of observations targetted a 
field around the W5-E Hn region. This region was se- 
lected because of its simple geometry and, at 2kpc (e.g. 
iKarr fc Martini I2003T I. is one of the nearest region s of trig- 
gered massive star formation (Mcg eath et al.|[2008l ). A single 



selected region. The available 850 fim array typically only 
had 50 per cent of the detectors working (although the num- 
ber used in the data reduction is selected dynamically and 
hence varies with time). Calibration was performed using an 
internal flat-field source, as well as absolute measurements 
using known calibrators on the sky. 



2 'Pong' is the default scanning mode for covering large areas 
with SCUBA-2, and is like a raster-scanning pattern except that 
it is designed to visit different parts of the map on different 
timescales. 'Rotating pong' means that the orientation of the pat- 
tern is allowed to rotate in sky coordinates as the observations 
are carried out, resulting in scans at many different angles, which 
is better from a map-making perspective. 
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region of 1° x 1° was mapped on 2009 December 5, at a speed 
of 600 arcsec s _1 (the nomi nal scan-rate for SA SSy) using 
the 'pong' scanning mode l|Kacklev et alJl2010h . The scan 
pattern had each successive sweep separated by 120 arcsec, 
ensuring an overlap to improve mapping performance. 

The on-source time was 70 minutes which resulted in 
a total of 12 passes over the entire region. In addition the 
central 0.5° x 0.5° was mapped at the same speed for an- 
other 70 minutes, covering that region 42 times, for a total 
observing time of 140 minutes. The average integration time 
per pixel for the inner and outer regions were 18 and 3 sec- 
onds respectively. The 225-GHz optical depth varied from 
0.08 to 0.12 with a mean of 0.10 (corresponding to Grade 3 
weather). This is substantially lower than the allotted opac- 
ity band for SASSy; however, the central 0.25 deg 2 region 
of the map has a sensitivity slightly exceeding the SASSy 
target level, and therefore represents a valid test of the de- 
tectability of sources. The central portion of the map has 
a noise of 25mJybeam _1 ; the value for the outer region is 
60mJybeam~ 1 , as reported by the noise map. Note that 
450 ;Um data were also obtained and were also reduced, but 
were found to be of limited use, with only the single bright- 
est source being detected. 
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Figure 3. Histogram of matched-filtered pixel values in the 'ex- 
tragalactic' map. This can be used to determine the spatial noise 
in the map. A 25.5 mJy noise curve is plotted for comparison, 
while the expected uncorrelated noise in the map is significantly 
smaller. 



3 DATA REDUCTION AND SOURCE 
DETECTION 

The raw timeseries data were process ed using smurf 
(|Chapin et al.l 120101; iJen ncss et alj l2010h called from the 



ORAC-dr pipeline |Gibb fc Jenness|[2010l '). smurf solves for 
the astronomical signal using an iterative technique, fitting 
and filtering out noise contributions from the atmosphere 
and the instrument. Details of the map-maker may be found 
in Chapin et al. (2010, in preparation) . Readings from the 
JCM T water- vapour monitor (WVM, iDempsev fe Friberd 
2008) were used to correct for atmospheric extinction. The 
ORAC-dr pipeline was used to mosaic the individual obser- 
vations using inverse- variance weighting. 

In the NGC 2559 raw map it is clear that time- 
dependent noise shows up as large-scale structure, obscuring 
the detection of the galaxy. This comes from a combination 
of residual sky fluctuations, as well as oscillations inherent 
to the instrument when these data were taken. Fig.[T]shows 
the NGC 2559 field after smoothing with the beam (to im- 
prove the signal-to- noise ratio for sources). The resulting 
map is dominated by structure at roughly the scale of the 
SCUBA-2 array. However, since the core of the galaxy shows 
up as a compact source, we are able to use a point source 
filter to make an unambiguous detection. The 'matched- 
filter' method implemented within ORAC-DR subtracts the 
map smoothed with a larger Gaussian (30 arcsec in this 
case) from the map convolved with a Gaussian equal to the 
JCMT beamsize (14 arcsec at 850 /im), thus giving a 'Mex- 
ican hat' type spatial filter. The resulting map is able to 
enhance sources that are approximately point-like (since we 
expect extra-galactic sources found by SASSy to be approx- 
imately point-like). We experimented with different choices 
for filter shape, and found that the default within ORAC-DR 
is close to the best we can do in terms of signal-to- noise ratio. 
In principle we could improve things by using specific knowl- 
edge of the expected shape of NGC 2559, but that would not 



be helpful in a blind SASSy search (which is what we are 
preparing here). 

Fig. fj shows the NGC 2559 field after applying the 
matched-filter algorithm to remove the large-scale struc- 
ture. NGC 2559 is now plainly visible in the centre of the 
map. Fig. [3] shows a histogram of the pixel values found 
in the matched-filtered map of NGC 2559. After removing 
the background variations in the map, the spatial noise is 
found to be approximately 25.5 mJy (equivalent for a point 
source). This is larger than the average noise of 19mJy cal- 
culated by the data-reduction pipeline (i.e. given by the 
noise map), due to not having filtered out all the spatially 
correlated noise. Since shallow extragalactic fields should be 
composed of white noise plus a few sources, the signal-to- 
noise matched-filtered map is re-normalized to have an rms 
of unity before searching for sources. In other words we are 
effectively using the spatial rms noise, rather than the purely 
white noise estimate which comes from the pipeline noise 
map. 

Fig. U shows the W5-E Galactic star-forming region 
map, while Fig.[S]shows the related signal-to-noise map and 
Fig. [5] shows an annotated version with several previously 
catalogued regions labelled. Since the purpose of SASSy is 
to find and catalogue new sources for future study, it is sim- 
plest to ignore source structure and to filter the map with 
the matched-filter. For objects within the Galactic strip, we 
use a 60 arcsecond background subtraction for the matched- 
filter. This scale is chosen since we expect Galactic sources to 
be more extended than extragalactic sources. Fig. 1 101 shows 
the W5-E map after being matched-filtered. The central 
~ 0.5° x 0.5° region has a spatial noise level of 15 mJy, while 
the outer region has a noise level of approximately 30mJy, 
once processed by the matched- filter (see Fig. [6j> . Note that 
these values are different from the noise values in the unfil- 
tered map since they only apply to the detection of compact 
sources (and so should not be considered as being a noise 
'per beam', as would be usual for extended structure). 
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Figure 4. The W5-E star forming region as mapped by SCUBA-2 at 850 £im, smoothed with a Gaussian with a FWHM of 14 arcsec. 
Colour scale ranges from —100 to +500 mjy beam -1 . The central roughly 1/4 of the map can be seen to have lower noise than the rest. 
Known objects are labelled, as is the approximate centre of the W5-E Ho region. 



4 DETECTIONS 
4.1 NGC 2559 Field 

In order to extract sources from the matched-filtered maps 
in a manner which could be easily autom ated for the full sur - 
vey, we use the fellwalker algorithm (|Berrv et al.ll2007f ) . 
implemented in the STARLINK software package CUPID. We 
use this algorithm simply to associate contiguous blocks of 
pixels with a single source - the brightness is estimated from 
the peak value in the match-filtered signal-to-noise map. 
Given the simple nature of the source detection, the choice 
of algorithm is not critical. A minimum number of pixels of 
7 and a low RMS level are used in fellwalker to recover 
as many real and noise peaks as possible for producing his- 
tograms. 

Fig. [7] shows a histogram of signal to noise peaks of 
all the sources within the NGC 2559 field. NGC 2559 is 



clearly detected at a signal to noise of about 7, and there 
are two other candidate sources at a signal to noise around 
4.5. There are no obvious counterparts for these two sources 
in any relevant survey. The histogram shows that the de- 
tection of NGC 2559 is an outlier in the distribution, but 
the two other candidates seem consistent with the noise dis- 
tribution, and in a map with > 10 4 beam-size pixels they 
are not very unlikely. Neverthele ss we checked f or cou nter- 
parts in the 20-cm NVSS survey ICondon et all (|l998l ). We 
additionally extracted archival 20-cm VLA data of approxi- 
mately the same depth (lower bandwidth but larger integra- 
tion time), re-reduced them and added them to the NVSS 
data. No radio sources are detected in the combined 20-cm 
image at the positions of the two peaks in the SASSy map, 
where the noise level is 320-350 /iJy. 

We also explored different choices for the filtering, and 
found that although the detection of NGC 2559 is robust, 
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Figure 5. Signal-to-noise ratio image for W5-E. The colour scale ranges from —2 to +10. The 'cold core' referred to in the text is circled. 
Known objects are labelled, as is the approximate centre of the W5-E Hll region. 



the next most significant peaks vary in position and bright- 
ness. The SASSy survey plan is to carry out short follow-up 
observations of such candidates to distinguish between real 
objects and false positives (whether just noise excursions or 
mapping artefects). This pilot survey suggests that the level 
at which follow-up will be worthwhile is around the 5<r level 
for large maps. 

Fig. [8] shows an optical image (from the Digitized Sky 
Survey) of NGC 2559 with matched-filtered SCUBA-2 con- 
tours overlaid. NGC 2559 shows up at a flux density of 
156 ±26 mjy and would be detectable in data representative 
of SASSy (with target rms of 30mJy). 

4.2 W5-E Field 

Turning to the Galactic pilot map, this field contains a 
number of known sources, including three bright-rimmed 
clouds (SF012, 13 and 14), and the Hll region Sh 2-201. 



Note that SFO 14 contains the massive young stellar ob- 
ject AFGL 4029. The three bright-rimmed clouds have all 
been detecte d previously at sub millimetre wavelengths us- 
ing SCUBA (|Morgan et alj|2008h . 

All of these sources are detected in the SCUBA-2 map 
with good to high significance. A total of 27 sources were 
identified at >6-<j after applying the matched-filter, most 
of which are unknown at submillimetre wavelengths. The 
central portion of the map containing the W5-E H n region is 
devoid of dust emission. Table[T]lists the objects found with 
a signal-to-noise ratio greater than 6 using the matched- 
filter method. Of the objects in Table[T] 11 are brighter than 
150 mjy and would be detected by a blind SASSy survey at 
more than 5<r. Fig. [9] shows a histogram of signal-to-noise 
peaks of all the sources extracted from the W5-E map. 

By inspecting the unfiltered map it is clear that some of 
the 'sources' found in this way are parts of extended filamen- 
tary structures. Fig|4] shows a number of extended, filamen- 
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Figure 6. Histogram of matched-filtered pixel values in the W5- 
E map for f 00 square pixel blank sections of the inner and outer 
regions of the map. Noise curves were fit and determined to be 
15 and 30mjy, respectively. 
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Figure 7. Histogram of peak signal-to-noise values within the 
matched-filtered NGC 2559 field. NGC 2559 itself shows up at 
around 6.6<r, while the 2 candidate sources at about 4.5o" appear 
to be noise bumps (note that a 4<r event is not very unlikely, given 
the number of pixels in this map). 



tary features of low signal-to-noise ratio. While faint, they 
are undoubtedly real as they show good agr eem ent with the 
CO da ta presented bv lKarr fc Martini (120031) andlNiwa et all 
l|2009h . as well as Spitzer MIPS images (|Koenig et al.ll200sl v 
However, in practice, mapping such low surface-brightness 
features is beyond the scope of SASSy and falls into the 
realm of followup observations triggered by detecting new 
compact sources. The most interesting Galactic sources 
found by SASSy will be relatively isolated, and they will 
be discovered through applying a simple automated source- 
extraction procedure similar to that used here. Nonetheless 
it is encouraging to find that SASSy is capable of detect- 
ing extended features and there are many approaches which 
can be taken to characterise such morphology. Further ex- 
perience with SASSy data will guide our approach. 
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Figure 8. SCUBA-2 850 (im 5 and 7<r contours overlaid on an 
optical 3-colour image of NGC 2559 (derived from DSS data). 
The black cross marks the position of the NVSS radio source. 
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Figure 9. Histogram of peak signal-to-noise values within 
the matched-filtered W5-E field. The brightest W5-E source 
(AFGL 4029) has a signal-to-noise ratio of over 100 and is not 
shown. 



Of all the new submillimetre detections perhaps the 
most striking source in the SCUBA-2 map is that labelled 
as the 'cold core' in Fig. [5] It shows up clearly in the signal- 
to-noise ratio map at a significance of 7<j. Comparison with 
Spitzer data reveals that while it is also detected at 160 and 
70 /im, it is completely absent at 24 ^m. As shown in Fig. 1111 
the SCUBA-2 source lies within the boundary of a bright- 
rimmed cloud, externally illuminated by the O-stars in the 
W5-E cluster. This source appears to be cold, a conclusion 
which is confirmed by our analysis below. 
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Figure 10. The matched-filtered map of W5-E with circles indi- 
cating the > 6cr sources listed in Table [T] 
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Right ascension 

Figure 11. SCUBA-2 contours overlaid on Spitzer 24 fim image. 
The source labelled 'Cold core' is circled and shows up clearly in 
the submm, yet within a dark region of the Spitzer image. It is, 
however, detected at 70 and 160 fim with Spitzer. 

5 PROPERTIES OF NGC 2559 
5.1 Ancillary data 

We complement the SASSy observations of NGC 2559 with 
available archival mid- and far-IR data. NCG2559 has been 
observed by both IRAS and Akari, providing spectral cov- 
erage from 9-160 fim. These data can be used to constrain 
the spectral energy distribution (SED) of the source. 

We estimate IRAS fluxes using the IRAS Scan Pro- 
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02 
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02 


57 
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03 
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33 
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Table 1. List of objects found with a peak signal- to- noise greater 
than 6 in the W5-E region using the matched-filter method. 

cessmg and Integration Tool (SCANPI0). Standard reduc- 
tion parameters were used, tak ing into account the observed 
size of the galaxy ( ~ 3 arcmin, IPrugniel fc Heraudeau|[l998l : 
IJarrett et ail [2003) as an additional constraint. Although 
NGC 2559 is actually detected as a point source in all IRAS 
bands, we assume a minimum distance of 6 arcmin from the 
nominal position of the source for background subtraction. 
We obtain Akari fluxes at 9 and 18 /tm fro m the ^4fcan/IRC 
All-S ky Survey Point Source Catalogue (jlshihara fc et all 
l2010h . and at 65, 90 and 140 ^m fro m the Akari/FIS All-Sk y 
Survey Bright Source Catalogue (jYamamura et aD l2010h . 
We decide not to use the 160 /im measurement owing to 
the extremely high noise affecting this band: this is con- 
firmed by the \ 2 increasing by about a factor 2 when this 
measurement is included in our fits. In addition to that, we 
estimate an upper limit on the 450 [im flux from the shorter 
wavelength SCUBA-2 map. A 20 per cent uncorrelated un- 
certainty is applied to all data. The observed mid- to far-IR 
data of NGC 2559 are shown in Table O 

5.2 SED fitting: dust models 

We fit the dust models of lDraine fell (|2007f ) to the com- 
pilation of available data. These models provide the dust 
emissivity per hydrogen atom, >(<?pah, t/ min , U m&K ), which 

3 Available at http://scanpiops.ipac.caltech.edu/9000/- 
applications / Scanpi / index, html 

4 Available at http://www.astro.princeton.edu/~draine/dust/- 
dust.html . 



A Pilot Study for the SCUBA-2 All-Sky' Survey 9 



Table 2. The observed SED of NGC 2559. Flux densities and 
respective errors are in Jy. Quoted errors include an additional 
20 per cent calibration uncertainty. 



A ((im) 
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SS U 


Instrument 
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1.5 


0.3 


Akari/IRC 


12 


1.4 


0.3 


IRAS 


18 


l.S 


0,1 


Akari/IRC 


25 


2.8 


0.6 


IRAS 


60 


2(5 


5 


IRAS 


65 


23 


5 


Akari/FIS 


90 


.39 


8 


Akari/FIS 


100 


6(5 


1.3 


IRAS 


140 


52 


12 


Akari/FIS 


450 


2.3 


2.3 


SCUBA-2 (u.l.) 


850 


0.156 


0.04 


SCUBA-2 



is a function of three parameters: the fraction of dust mass 
in PAHs, qpAW, the intensity of the radiation field from stars 
heating the interstellar medium, I7 m in; and the intensity 
of the radiation field in photodissociation regions (PDRs), 

Umax ■ 

We apply th e same general method explained in 
iDraine fe Lil (|2007m , b ut we use on l y the seven Milky Way 
dust model sets, as in lWiebe et al.l l|2009l ). Accordingly, we 
also set Umax = 10 6 (we tested the reliability of this as- 
sumption by leaving Umax as a free parameter and found 
that the fit returns the same value). We thus fit to the ob- 
served SED a linear combination of diffuse ISM models (with 
t/max = f/min) and PDR models: 

-Fi/(gPAH, Umin, Umax) OC 

X K 1 _ 7)>(<?PAH, Umin, Umin) 
+ 7>(<?PAH, Umin, Umax)] , 

where Md is the dust mass, jtih is the mass of a hydrogen 
atom and D = 20.8 Mpc is the distance to the galaxy. We 
then use the derived values to evaluate the dust-weighted 
starlight intensity ((/). The best-fit values for the parame- 
ters are found through \ 2 minimization. Fig. 15.21 shows the 
observed SED and best-fit curves. 

We notice that the Akari/FIS data are systematically 
low with respect to the IRAS data. We thus evaluate the 
outcome of the fit using either the full data set or a subset 
to assess the effect of this offset. Table Rx2"l summarizes the 
different fits and results. 

Fitting the dust models to the full set of data points 
available yields a reasonable fit, although the \ 2 is relatively 
high. We obtain the values: (/pah = 3.2; 7 = 0.01; Umin = 20; 
and (U) = 21.9. We derive a dust mass of 1.5 x 10 7 M Q and a 
far-IR luminosity of 2.0 x 10 10 Lq, from which we calcul ate a 
star-fo rmation rate of 2.8 Mq yr _1 using the relation of lBelll 
(2003). The derived dust mass is about 2 times larger than 
the value obtained bv lBettoni et all |2003l ), who find a value 
of 8.3 x 10 6 M Q from data at 60 and 100 /im. 

Note that the SFR can also be estimated using the 
20 cm radi o flux density along with the FIR/radio corre- 
lation (e.g. ICondon|[T992l : ICram et al.lll998l : iHopkins et all 
120011 ). With 51.4 = 260 mjy this gives an SFR of around 
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Figure 12. The observed far-IR and sub-mm SED of NGC 2559, 
together with dust model fits. An arrow marks the 3cr upper limit 
at 450 /tin. The solid (black) line shows the best fit to all available 
data points. The dotted (blue) line is the best fit to the datasct 
excluding the Akari/FIS data. The thin solid (black) line is the 
best fit to the dataset excluding the IRAS data. The dashed 
(green) line is the best fit model without the SCUBA-2 points. 
The dot-dashed (orange) line is the modified blackbody fit with 
/3 = 2. The triple-dot-dashed (red) line is the modified blackbody 
fit with free /3. Error bars include an uncorrelated 20 per cent 
uncertainty. 

lOMQyr -1 , which is fairly consistent with the results ob- 
tained from the FIR fit. 

Exclusion of the Akari/FIS data points yields again a fit 
which correctly matches the SCUBA-2 850 /mi point. Best- 
fit parameters are rather similar, the only noticeable change 
being the increase of t/ m i n = 25 and of (U) = 27 A, with the 
resulting dust mass, FIR luminosity and SFR fairly consis- 
tent with the previous fit. 

Removing the IRAS points from the data set yields a 
similar fit with a \ 2 of 10.7 with 3 degrees of freedom. The 
fit parameters are almost undistinguishable from the first fit 
at A ^ 30 /im. 

For comparison, we also fit the models after excluding 
the SCUBA-2 point and upper limit. We obtain lower values 
of t/min = 5 and (U) — 5.6). But now the best-fit misses the 
850 /im flux density by about a factor 4. 

This shows the strong leverage of the SCUBA-2 data to 
properly constrain the FIR and sub-mm SEDs of galaxies. 
In particular, in addition to the dust mass, the value of the 
starlight intensity £/ m in (and consequently {U}) is strongly 
sensitive to the 850 /im flux density. 

5.3 SED fitting: modified blackbody 

For comparison with the more detailed IDraine fc Li (|2007l ) 
models, we also fit a modified blackbody spectrum to the 
observed data. The fit is carried out assuming a shape for 
the modified blackbody described by the expression 

where vq = 1.2 THz = c/(250/tm), k is the dust mass ab- 
sorption coefficient at vq, /3 is the dust emissivity index and 
Md is the dust mass. Once again, the best-fit parameters 
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Table 3. lDraine fe~Ll l|2007l l model parameters for NGC 2559 and derived physical quantities. 



Data removed 


X 2 


degs. of 


<?PAH 


^min 


(U) 


7 


M d 




SFR 






freedom 










(1O 7 M ) 


(1O 1O L ) 


(Moyr" 1 ) 


None 


17.7 


7 


3.2 


20 


21.9 


0.01 


1.5 


2.0 


2.8 


Akari/FIS 


10.8 


4 


3.2 


25 


27.4 


0.01 


1.3 


2.1 


3.0 


IRAS 


10.7 


3 


4.6 


20 


22.0 


0.01 


1.5 


1.9 


2.7 


SCUBA-2 


3.5 


5 


4.6 


5 


5.6 


0.01 


5.7 


2.2 


3.1 



Table 4. Modified blackbody fits for NGC 2559 and derived physical quantities. 



Fit 


x 2 


degs. of 


T d 


M d 


£fir 


SFR 






freedom 


(K) 


(10 7 M ) 


(10 10 L o ) 


(Moyr- 1 ) 


Fixed /3 


6.0 


5 2 


29 


0.31 


1.8 


2.6 


Free 


4.7 


4 2.3 


26 


0.47 


2.0 


2.8 



are found by x 2 minimization. We a ssume a mean value of 
k — 0.29 (see e.g. IWiebe et al.ll2009l . although there is con- 
siderably uncertainty in this value) and fix f3 — 2. 

This modified blackbody fit yields consistent values of 
Md, £fir and SFR with respect to the detailed dust mod- 
els. Table 15.21 summarizes the modified blackbody results, 
and the fitted curves are shown in Fig. 15.21 We see that the 
dust temperature of NGC 2559 is around 26-29 K. This is 
warmer than for most of the galaxies studied in the SINGS 
sample iDraine fc et alj 1 20071). as well as those studied us- 
ing BLAST lWiebe et all (I2009T I. This is consistent with re- 
quiring a relatively higher value of U than for those other 
galaxies. 



6 SOURCE PROPERTIES IN W5-E 

From our catalogue of compact sources within the W5- 
E region we select two example sources, representative of 
two extreme regimes: (i) the brigh test source in the field, 
AFGL4029 l|Deharveng et al.|[l997l and references therein); 
and (ii) a fainter source to the NW of the first, as an example 
of a potentially colder core (labelled 'cold core' in Fig. [5}. 
The SCUBA-2 flux densities for both sources were obtained 
from photometry within an aperture of 30-arcsec diameter. 

Ancillary data available for this region include Spitzer 
MIPS 70 /im observations of the whole field, as well as IRAS, 
Akari coverage, and SCUBA imaging for part of the field. 
The SCUBA-2 data are in excellent agreement with the 
SCUBA measurements, being consistent within errors, but 
covering a much wider area. We show Spitzer 70/im contours 
on top of the SCUBA-2 image in Fig. [13] 

The bright source AFGL4029 is detected in almost all 
wavebands, having upper limits only at 65 and 90 fim. We fit 
a modified blackbody spectrum to the observed SED, after 
adding a 20 per cent uncorrelated uncertainty to the errors, 
as we did for NGC 2559. Leaving /3 as a free parameter in the 
fit yields a value of /3 = 2.02, which suggests that /3 = 2 is a 
good assumption for the dust emissivity index. We derive a 
temperature T = 31 K and a dust mass of 60 Mq, with a FIR 
luminosity of 5000 Lq. This is in r easonable agreement with 
the value of 3200 L estimated bv lMorgan et all (|2008l ). 

The faint source is detected only at 70 and 850 /im. A 



relatively nearby 90 ^m source is detected by Akari, however 
its position is offset by 28 arcsec from the Spitzer position. 
Although this lies within the PSF of Akari at 90 /im, we 
prefer not to use this measurement in our analysis to avoid 
contamination by other potential sources (and in any case it 
adds little to the 70/im constraint). Reliable flux densities 
were obtained from the Spitzer 24 and 70 images, using 
the same 30-arcsec aperture as for the SCUBA-2 data. While 
the cold core is also evident in the 160 /im MIPS image, the 
presence of significant artefacts prevents the extraction of a 
useful flux measurement. 

Although only two photometry points are available, a 
fit to the SED using again a modified blackbody with /3 = 2 
yields reasonable quantities and confirms this source to be a 
cold protostellar core. We obtain a temperature T = 17 K, 
FIR luminosity of 27 L and a dust mass of 11 M . Com- 
paring the Akari 90 fim point to the fit shows that the mea- 
surement is in excellent agreement with the expected flux 
density, suggesting that this might actually be the counter- 
part to the same source. The result of the fit for both sources 
is shown in Fig. [S] 

The readiness with which we were able to pick out a rel- 
atively cold source, even in this small pilot study, shows that 
SASSy will be able to detect single low-mass cold clumps in- 
side larger star-forming regions. Detailed follow-up of such 
sources will determine where they lie on the star-formation 
sequence. 



7 DISCUSSION 

With the current SASSy S2SRO observations, we have 
shown that we will be able to detect and catalogue real 
sources, both Galactic and extragalactic. Based on early 
commissioning estimates, SCUBA-2 mapping during S2SRO 
was expected to be approximately 50 times faster than for 
SCUBA. A direct comparison with the SCUBA data in the 
W5-E field shows the mapping speed improvement factor 
to be 60, after scaling the map sizes and noise levels. As- 
suming similar weather conditions and a similar efficiency 
for the new science-grade arrays currently being commis- 
sioned on SCUBA-2, SASSy should be able to map the sky 
at a rate of > 0.8 deg 2 per hour to the target sensitivity, 




Right ascension 

Figure 13. SCUBA-2 850 /im emission in colour with Spitzer 70 jim contours overlaid. Sources are labelled as in Fig. [4] The Spitzer 
contours do not close to the south of W5-NE due to a lack of coverage there. 




100 

Wavelength [^im] 



100 

Wavelength [^m] 



Figure 14. The observed far-IR and sub-mm SED of AFGL 4029 and of the cold source in W5-E. Black points are used for the fit, 
arrows mark upper limits. The dotted (red) line is the modified blackbody fit with f) = 2. 
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exceeding the capabilities of SCUBA by hundreds of times. 
However, the performance of the new 850 fj,m arrays are cur- 
rently unknown, since they are yet to be tested on the sky, 
and hence the area which will be ultimately mappable by 
SASSy is still uncertain. At the present time SCUBA-2 has 
all eight arrays installed and is undergoing the first stages of 
full commissioning. It is anticipated that this will be com- 
pleted in mid-2011 and the instrument made available to 
the community as soon as possible thereafter. SASSy will 
therefore start in earnest some time in 2011. 

Nevertheless, these initial data have shown that it 
is possible to reach a la sensitivity to point sources of 
~ 30 mjy while scanning rapidly with SCUBA-2 in relatively 
mediocre weather conditions. The map-making and source 
extraction procedures are fast and require little in the way of 
human intervention. Experience with this pilot programme 
has already been fed back directly into the software pipeline 
and most of the data processing and analysis is now auto- 
mated. 

The science case for SASSy remains strong and has been 
made stronger by recent Herschel discoveries. There is still a 
pressing need for a wide-area shallow 850 /im survey in the 
era of Herschel, Planck and ALMA. This pilot study has 
shown that it is feasible to find sources in the shallow maps 
which SCUBA-2 will soon produce routinely. 
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